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The  subject  of  this  paper  is  a  version  of  the  inverse  reflection 
problem  for  a  smoothly  stratified  elastic  medium.  The  same  mathematics 
describes  the  inverse  problem  of  the  vibrating  string.  This  problem 
is  solved  in  a  constructive  way.  Also,  a  priori  estimates  are  de¬ 
rived  which  exhibit  the  continuous  dependence  of  the  solution  (index 
of  refraction,  relative  sound  speed)  on  the  data  (scattering  or  reflec- 
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SIGNIFICANCE  AND  EXPLANATION 


In  an  inverse  reflection  or  scattering  problem,  the  mechanical 
(say,  elastic)  properties  of  a  medium  are  to  be  determined  from  knowl 
edge  of  incident  and  reflected  waves.  Such  problems  arise  in  explor' 
ation  geophysics,  materials  testing,  plasma  diagnostics,  physical 
chemistry,  ultrasound  radiology,  and  many  other  contexts.  A  simpli¬ 
fied  problem  of  great  practical  importance  is  the  inverse  acoustic 
reflection  problem  for  a  smoothly  stratified  medium  in  which  the 
density  varies  as  a  function  of  depth  only. 

This  report  gives  a  constructive,  exact  solution  of  this  problem 


using  only  techniques  available  also  for  the  three-dimensional  inverse 


acoustic  problem  for  unstratified  elastic  media.  Additionally,  the 


methods  used  here  provide  stability  estimates  which  are  extremely 


important  for  applications  and  for  numerical  implementation 


The  responsibility  for  the  wording  and  views  expressed  in  this 
criptive  summary  lies  with  MRC,  and  not  with  the  author  of...th 
report. 


THE  INVERSE  REFLECTION  PROBLEM  FOR  A 
SMOOTHLY  STRATIFIED  ELASTIC  MEDIUM 

W.  Symes 

§  1 .  Introduction 

The  subject  of  this  paper  is  a  version  of  the  inverse  reflection 
problem  for  a  smoothly  stratified  elastic  medium.  The  same  mathematics 
describes  the  inverse  problem  of  the  vibrating  string.  This  problem  is 
solved  in  a  constructive  way.  Also,  a  priori  estimates  are  derived  which 
exhibit  the  continuous  dependence  of  the  solution  (index  of  refraction, 
relative  sound  speed)  on  the  data  (scattering  or  reflection  measurements). 

We  talte  particular  care  to  use  only  constructions  which  apply  in 
principle  to  higher -dimensional  problems  of  a  similar  sort.  Most  other 
treatments  of  this  problem  proceed  via  a  reduction  to  an  inverse  Sturm- 
Liouville  problem  (see  e.g.  [1]  and  (2)).  This  coordinate  transformation 
is  no  longer  available  in  higher-dimensional  problems,  so  we  avoid  it, 
except  where  it  can  be  interpreted  as  a  coordinate  transformation  along 
the  rays  of  geometric  optics;  this  is  the  case  for  the  a  priori  estimates 
of  c  at  the  end  of  Section  4.  (See  ?4  for  a  discussion  of  this  point). 

Other  authors  who  avoid  dependence  on  peculiarly  one-dimensional 
tricks  have  invariably  used  approximate  methods  (JWKB,  Born  series:  see 
[3],  (4),  Ch.  XIII,  and  practically  any  article  in  the  literature  of  inverse 
scattering  in  exploration  geophysics,  physical  chemistry,  ultrasound 
radiology,  etc.)  In  contrast,  our  methods  are  exact:  we  construct  the 
index  of  refraction  exactly  by  an  iterative  approximation  procedure. 

Our  formulation  of  the  inverse  reflection  problem  is  time-dependent, 
which  also  contrasts  with  many  other  treatments  (e.g.  [2],  [5]).  The 
solution  propounded  here  can  also  be  adapted  to  steady-state  (frequency- 
domain)  problem  formulations,  other  boundary  conditions,  etc. 

We  abjure  further  discussion  of  the  literature  and  proceed  to  describe 
our  results. 

Sponsored  by  United  States  Army  under  Contract  No.  DAAr,29-75-C-0024 
This  material  is  based  on  work  supported  by  the  National  Science 
Foundation  under  Grant  No.  MCS78-09525. 
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The  physical  setting  of  our  problem  is  as  follows.  Consider  a 
stratified  elastic  medium  with  unit  elastic  moduli  (L.ame'  constants). 

(This  latter  restriction  may  be  removed,  in  which  case  we  recover  the 
elastic  impedence  rather  than  the  acoustic  index  of  refraction).  The 
density  p  varies  as  a  function  of  only  one  of  the  coordinates,  say  x. 

We  assume  that  p  is  known  and  constant  in  the  half-space  (x  ^  O)*  In 

fact,  for  the  sake  of  simplicity  p  •  1  for  x  ^  0,  and  smooth  (of 

class  )  otherwise.  The  methods  can  be  modified  to  cope  with  dis¬ 
continuous  p  —  that  is  easy.  The  hard  part  is  smooth  p  .  We  set  up 
a  measuring  device  at  x  =  O,  and  record  the  (infinitesimal)  reflected 
wave  F(t),  t  >  O,  resulting  from  an  (infinitesimal)  impulsive  incident 
wave  of  the  form  6(x-t),  t  <  0,  impinging  on  the  unknown  medium 
lx  0l. 

Our  results  are  phrased  in  terms  of  the  index  of  refraction 
-1/2 

c  =  p  .We  obtain  sufficient  conditions  on  F  (which  are  very  close 
to  necessary  —  see  [12]),  as  in: 

Theorem  1.  Suppose  T  >  0,  F  :  [0,2T)  H  i£  class  C^, 

F(0)  c:  O,  and  the  kernel 

H(s,t)  E  iF(s  +  t)-f®  dT  F(8-T)F(t-T) 

2  •'o 

^  ^  9  <:  t  ^  T, 

H(t,s)  s  H(s,t) 

defines  a  self-adjoint  Hilbert-Schmidt  operator  M  on  L^[0,T]  with 
the  property 

H  +4H  N  e  >  O 

Then  there  exists  a  unique  x  >  0  and  a  unique  positive  function 
c  :  [0,y)  (c  '  0)  of  class  so  that  f  c“^  «  T  and  F(t)  <= 

- - •'o 

u(0,t),  t  ,  0  for  the  solution  u  of  the  initial  value  problem  (for 
which  define  cal  for  x  ^  0): 


Then  the 


whose  existence  is  implicit  in  the  above  state 


ment  is  continuous  (even  locally  Lipschitz)  as  a  map 


The  proof  proceeds  along  the  lines  laid  down  in  (6) ,  where  a  similar 


problem  was  solved  with  the  differential  equation  (1.1a)  replaced  by 


In  reference  [6],  therefore,  the  characteristics  of  the  problem  are  known 


from  the  outset.  The  characteristics  are  exactly  what  is  to  be  found 


in  the  present  problem,  however,  which  gives  it  the  nature  of  a  free 


boundary  problem.  This  nature  is  unavoidable  in  higher-dimensional 


inverse  problems  for  elastic  waves,  and  we  meet  in  head-on,  which  causes 


more  difficult  to  derive  than  are  the  corresponding  estimates  in  [6J 


The  main  tool  is  the  progressing  wave  decomposition  of  the  solution 
u  of  an  initial  value  problem  (1.1): 


t 


i«  the  travel-time  function,  and  K  is  smooth  inside  the  light  cone  with 
apex  (0,0)  (Section  2). 

This  expansion  is  established  in  §2  under  a  condition  on  c  which 

guarantees  that  distribution  solutions  of  singular  initial  value  problems 

2 

i  ^  2  A 

for  the  operator  □  »  -  c  — ^  are  unique.  This  condition  is 

dt^  ax'' 

not  necessary.  It  is  proved  in  ref.  (12)  that  a  sharp  condition 
is  c  €  .  Also,  a  converse  to  Theorem  \  is  given  in  (12),  which 

amounts  to  a  complete  solution  of  the  inverse  problem  (necessary  and 
sufficient  conditions  on  F). 

In  §3,  it  is  shown  that  the  pair  (K,c)  solves  the  GL  system  of 
Volterra  equations  (Theorem  B) : 

H(T(x),t)  «.  2  c~^  '^(x)Kvx,t) 

+  dyc"^t,y)K(y  ,T(x)  )K(y  ,t) 

O 

for  t  N  T(x)  (1.2a' 

K(T(x),x)  -  ^  (1.2b) 

The  initials  GL  stand  for  Gel'fand  and  Levitan,  for  they  introduced  inte¬ 
gral  equations  of  this  sort  into  inverse  scattering  theory  in  their 
seminal  paper  (13).  Our  wor)c  is  in  direct  line  of  descent  from  theirs: 
the  "nonlinear  integral  equation"  of  (13)  is  derived  by  the  present 
techniques  in  (6) ,  GL  also  stands  for  q roup  law t  indeed,  the  equation 
(1.2a)  expresses  in  compact  form  the  group  property  of  the  solution  oper¬ 
ators  for  the  Cauchy  problem  for  which  follows  from  the  time- 

independence  of  its  coefficient. 

A  number  of  crucial  a  priori  estimates  are  determined  in  Section  4 
for  the  solution  (K,c)  of  (1.2).  These  involve  sup  norms  of  F  and  the 


5 


number  This  latter  number,  although  rn  principle  present  in  the 

data  F,  is  in  practice  difficult  to  extract.  On  the  other  hand,  in 
practical  problems  one  often  has  known  bounds  on  the  index  of  refraction. 

It  is  shown  in  Section  6  that  a  priori  bounds  on  c  and  its  first  two 
deiivatives  determine  a  lowei  bound  for  e  ,  hence  can  be  employed  in 
place  of  6  in  the  a  priori  and  stability  estimates.  We  use  the  results 
of  Section  4  in  Section  5  to  show  that  an  iteration  scheme  converges  to 
a  global  solution  of  the  GL  system  (1.2).  The  solution  defines  a  continuous 
map 

7.J,  -  r  =  {{X.c,K)  ;  X  >  O}, 
c  t  C^{0,xj  ,  K  fc  C^(2.(T.c))  I 

where 

(J(T,c)  »  {(x.t)  X  ^  X.  Tlx)  ^  t  ^  2T-T(x)1. 

Finally,  we  establish  in  Section  7  that  the  solution  (K,c)  of  the  GL 

system  actually  solves  the  Chudov  laoundary  value  problem 
2  2 

(-Ay  -  c^  -A)K  =  O  in  ^(T,c)  (1.3a' 

dx‘‘ 

K{0,t)  =  F(t),  (^  -  ^)K(0,t)  s  0  (1.3b' 

V,  1  1/2,  1/2,  1/2.,  ,  ,  , 

K(x,T(x))  =  j  c  '  (x)J  c  (c  )  .  a. -'O' 

Since  it  was  established  in  Section  2  that  the  (necessarily  unique) 
solution  of  (1.3)  is  the  smooth  part  of  the  solution  of  the  singular 
initial  value  problem  (1.1),  it  follows  that  c  is  the  solution  of  the 
inverse  problem  sought  in  Theorem  A.  which  completes  its  proof. 

The  result  of  Section  7  also  establishes  that  the  GL  system  (1.2) 
and  the  Chudov  problem  (1,3)  are  completely  equivalent.  In  particular, 
the  a  priori  estimates  of  Section  4  also  hold  for  the  solution  of  (1.3'. 
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In  the  context  of  the  inverse  problem  of  16),  a  similar  observation  was 
used  in  17)  to  prove  stability  of  an  optimally  efficient  numerical 
scheme  based  on  a  Chudov  problem.  In  fact,  the  Chudov  problem  of  [6), 

(7)  was  suggested  as  an  approach  to  the  inverse  spectral  problem  for  the 

H 

Schrodinqer  operator  by  Chudov  in  the  mid  50's  (see  (14),  appendix), 
hence  the  name.  The  present  results  will  likewise  provide  the  base 
for  a  stability  result  and  consequent  a  priori  error  estimates  for 
efficient  numerical  solution  of  the  present  problem.  This  matter  will 
be  discussed  elsewhere. 

To  end  this  introduction,  we  note  that  the  present  results  can  be 
combined  with  well-known  techniques  from  exploration  seismology  to  solve 
the  inverse  problem  for  piecewise  index  of  refraction  with  jump 

discontinuities.  Finally,  we  mention  the  work  of  O.  Hald  [8]  on  inverse 
Sturm-Liouville  problems  and  P,  Deift  and  E.  Trubowitz  (9)  on  inverse 
Schrodinger  Scattering.  These  authors  also  prove  stability  results,  and 
their  techniques  are  actually  closely  related  to  those  presented  here. 


2  .  The  progressing  wave  expansion 

The  goal  of  this  section  is  to  express  the  solution 
singular  initial  value  problem 


-  u  »  -  c^ 


r)u 


ax' 


U(X.O)  =  i  (X  -Xq) 


u  of  the 


(2.1a) 

(2.1b) 


—  (X.O)  s  O 

at 


as  a  superposition  of  progressing  waves 

Ug(x,t)  =  2  c^^^{x)c"^''^(Xq) 16 (t  +  T(x.Xq) )  +  6 (t  -  T(x.Xp) ))  (2.2) 

where 

X 

T(x,x.)  =  /  c 

^0 


plus  a  remainder  K(x,t;XQ)  which  is  smooth  inside  the  light  cone  with 
apex  vx_,0)  and  has  a  jump  discontinuity  on  the  boundary  of  the  light 


■one  given  by 


The  first  step  is  the  determination  of  conditions  on  a  compound 
distribution  (see  [10],  Ch.  7) 


which  ensuies  that  C  v  has  singularities  no  worse  than  those  of  v 


in  this  case  that  no  derivatives  of  the  delta  distribution  appear 
y  along  the  surface  =  0.  Here  g  and  CP  are  understood  to  be 


smooth  functions,  with 


playing  a  role  analogous  to  that  of  the  phase 


function  in  geometric  optics 


A  short  calculation,  using  the  chain  rule  for  compound  distributions 


1)  the  eikonal  equation 


Two  convenient  independent  solutions  to  the  eikonal  equation  are 


passing  through  (x_,0) 


The  curves 


The  corresponding  transport  coefficients  must  solve 


to  first  order.  Now  demand  that  for  suitable  choice  of 


the  superposition 


should  satisfy  the  initial  conditions  (2.1b).  We  obtain 


The  last  is  in  fact  required  to  hold  to  first  order  at  x 


easily  checks  that,  if 


then  the  equations  (2.5)  are  satisfied,  including  the  last  one  to  first 
order,  provided  that  the  transport  equations  (2.4)  are  taken  into  account 
To  solve  the  transport  equations ,  notice  that 


must  solve 


along 


with  initial  condition 


The  initial  condition  is  satisfied  if  we  take 


r 


I 


order  iff  it  is  satisfied  identically  in  x  and  t,  in  which  case  the 
solution  is 


J 
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Finally,  one  shows  that  the  problem 


_  w  =  -Z  (u  +  u  ) 
— c  c '  s  c  ^ 


W  =  H  0  for  t 
0  t- 


O 


has  a  unique  continuous  solution  w  supported  in  the  light  cone  C(Xq). 
This  is  spelled  out  in  ([11],  Ch.  VI,  §4)  and  ([10],  Ch .  7)  under 
stringent  differentiability  assumptions  on  c,  and  in  ([12],  §2)  under 
weaker  restrictions  more  suited  to  the  discussion  of  the  inverse  pro¬ 
blem,  In  any  case,  we  can  now  represent  the  solution  u  of  (2.1)  as 


u  =  u  +  u  +  w. 
s  c 

Since  w  vanishes  on  3C(Xq)  and  is  continuous,  the  summand 
K  =  u^ + w  has  the  asserted  jump  along  3C(Xq) ,  by  virtue  of  (2.8). 


§  3 .  The  GL  equation 

Since  any  function  may  be  written  as  a  superposition  of  delta 
functions,  we  can  write  the  general  solution  u(x,t)  of  the  partial 
differential  equation 

□  u  =  O 
c 

satisfying  the  initial  condition 

l^(x.O)  .  0 

as 

00 

u(x,t)  =  J  dx^  S(x,t;x^^,0)u(x^,0) 

—CO 

where 

S(x,t;xQ,0)  =  u(x,t) 


is  the  solution  examined  in  the  previous  section:  that  is 


o 
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-c‘^'  ■ 


S  Ix.0;Xq,0)  «  Kx  -  Xq) 


D2S (x,0;Xq,0)  »  0. 

S  is  i elated  to  the  Riemann  function  ^  of  C!  .  which  solves  the 

-  g 

initial  value  problem 

“c'^t  '  ° 

R(x.0;Xq.0)  p  0 

D2R(x,0jXq,0)  -  6(x-Xq) 


(see  ( 11) ,  Ch.  V,  §5 )  by 


Define 


R  ( X , t ; Xq , 0 ) 


S(x,s;Xq,0) . 


R(x,t;XQ,tQ)  -  R(x,t  -  t^jx^.O) 

S  ( X ,  t  f  Xq  t  tQ )  *  S(x*t“  tQ  ?  Xq  #0)  (3*1) 

Then  the  general  solution  of  C^u  -  O  with  arbitrary  initial  data  is 
given  by 

u(x,t)  -  y  dXQlS(x,t;XQ.tQ)u(XQ,tQ)  +  RCx.tfXQ.t^)  ^  u(XQ,tQ)) 

— ®  0 


This  works  because  the  coefficiant  of  is  independent  of  t. 

For  the  same  reason,  the  initial  value  vector  (u(’,t  )j  is 

propagated  in  t  by  a  group  of  operators  (bounded,  in  fact,  in  a 
suitable  function  space).  The  distribution  kernel  of  this  group  of 
solution  operators  is  the  matrix 

S  R  \ 

D^S  DjR  ) 
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Tlie  group  law  is  expressed  in  terms  of  R  by 


<?(x,tfx^,tjj)  »  dy/^lx.tjy  ,t' 


(3.2) 


This  relation  is  fundamental  for  what  follows.  Since  the  compi'nents 
of  may  all  be  expressed  in  terms  of  the  scalar  kernel  S,  it  is 

plausible  that  (2.2)  may  be  expressed  in  terms  of  a  property  of  S. 
This  is  indeed  the  case. 

To  derive  this  relation,  we  require  some  further  symmetries 
of  S,  which  are  as  follows: 


1)  S  is  even  in  is  odd  in 


2)  The  kernel 


S  (x.tfXQ.t^)  =  S ( Xq, tpjx , t ) 


is  a  solution  of  the  adjoint  equation 

O 

(see  [llj  ,  Ch.  V,  §5.3) . 


*  * 

D  S 
c 


Now 


at^  dx^ 


So  (We  suppress  for  the  moment  the  dependence  on 

.2  .2 


O  »  ~  c^(x))S  (x,t)  ”  — ^  S  (x,t)  -  -^^(c^txlS  (x,t)) 

at  ax''  at  ax' 


(c"^(x) 


2 

2  -  -^■)(c^(x)S  (x.t))  =  c“^(x)Z  (c^(x)S  (X,t)). 

at  dx"^ 


2  * 

Since  c  >  0,  it  follows  that  c  {x)S  ix,t)  solves 


lZ^  c^(x)S  (x,t)  =  0. 


Now 


c^(x)S*(x,tQ)  *  c^(x)6(x-Xq)  •=  c^(x^,)a(x  -  Xj^) 


if 
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and 


A  -  0. 

We  assume,  without  further  explanation,  that  the  coefficient  c  is 
sufficiently  well-behaved  that  distribution  solutions  of  the  initial 
value  problem  for  the  wave  equation  are  unique  (conditions  on  c  under 
which  this  is  so  are  discussed  in  [12];  §2,  certainly  sufficient  is 
c  6  C^).  Therefore 


i.e. , 


S(Xo.to;x,t)  =  c“^(XQ)c^(x)S(x,t;XQ,tQ) 

which  is  the  desired  symmetry  relation. 

Now  examine  the  (1,1)  component  of  the  group  law  equation  (3.3) 
which  reads 

00 

S(x,trXQ,tQ)  =  J  dy  {S (x , t;y , t ' )S (y , t 'x^, t^) 

+  R(x,t;y,t')D2S(y,t';XQ,tj^)  K 

We  set  X  =  Xq,  =  O  and  replace  t  by  t  +  s,  t'  by  s  to 
obtain 


S(XQ,t+ s?Xq,0)  «=  J  dy  (S(XQ,t+ 8;y,8)S(y,s,XQ,0) 

+  R(XQ,t + s;y ,s )D2S (y,S ;Xq,0) ] . 

Using  the  time-translation  symmetry  (2.1)  rewrite  the  r.h.s.  as 
“  J  dy{S(XQ,t;y,0)S(y,sjXQ,0)  +  R(Xq, t;y, 0)028 (y,srXQ,0) 1 . 


Now  note  that,  since  S(y,8;xQ,0)  is  even  in  s,  its  s-derivative 
D2S(y ,srXQ,0)  is  odd.  The  second  term  in  the  integrand  is  therefore 
odd  in  s,  whereas  the  first  is  even.  Replace  s  by  -s,  add,  and 
divide  by  two  to  obtain 


dy  S(x  ,t;y,0)S(y,8jx_,0) 


Now  use  the  adjoint  symmetry  (2)  above  to  interchange  the  arguments 
in  the  first  factor  in  the  integrand: 


dy  c"  (y)S(y,0;x_,t)S(y,s;x  ,0) 


dye  (y)S(y,-t;x  ,0)S(y,s;x  ,0) 


dy  c~  (y)S^y,t;x_,0)S(y,s;Xf.,0) 


Conversely,  one  can  show  that  (3.3)  entails  the  group  law  equation  for 


the  full  Riemann  matrix 


One  can  show  ([12],  §2)  that,  under  suitable  smoothness 


to  obtain  a  continuous,  even  function  of  all  t.  Define  also 


S(x_,t;x_,0)  =  (5(t)  +  F(t) 


Recall  also  from 


(x_)  (6  (t  +  T(x,Xf^) )  +  6  (t  -  T(x,x_) )  ] 
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Since  Xp  will  remain  fixed  for  the  rest  of  this  discussion,  we  shall 
henceforth  set  Xq  =  0.  The  nuirtoer  c(Xq)  •  c(0)  is  now  surely  a 
positive  constant  which  can  be  set  equal  to  1  by  scaling  t,  which  we 
assume  has  been  done.  We  write 

T(x)  -  T(x,0)  «  c"^. 

‘O 

We  now  have  in  place  of  (3.3) 

2[4(8+t)+Ms-t)J  +  G(8,t)>=j  dyc“^(y)S(y,t}0,0)S(y,sj0,0)  (3.4) 

with 

S(y,t;0,0)  =  5  c^^^(x)  (6(t+ T(x))  +  6(T -T(x))]  +  K(x,t)  (3.5) 

One  now  substitutes  (3.5)  into  (3.4),  and  after  some  computation  eliminates 
the  singular  terms  to  obtain 

G(s.t)  =  ^(c"^^^(x"(s))K(x~(s),t)  +  c"^’^^(i<‘*'(s))K()c'^(3),t)J 

+  dy  c"^(y)K(y,t)K(y,s)  (3.6) 

'  X"(s  ) 

+ 

Here  s  «  Jf~(s)  «  x  is  the  inverse  function  to  xh  +  T(x)»s,  and 
is  thus  the  solution  of 

^  X"  =  +  c(X-),  X~(0)  =  O. 

The  light  cone  through  (0,0)  is  thus  described  by  ((x,t)  :  X~  ( t )  ^  x  ^x"*^  ( t ) ) , 
+ 

and  t  (X~(t),t),  X  (x,  +T(x))  are  equivalent  descriptions  of  the 
characteristic  curves  emanating  from  (0,0).  Recall  that  (Section  2,  (2.3)) 

K(x,T(x))  =  i  c^'^^(x)J'’‘  c^'^^(c^'^)'  x  N  o 
^  0 

K(x,-T(x))  -  T  c^/^(x)j’‘  c^/^(c^^^)'  x  i  O  (3.7) 

^  O 

In  view  of  the  formulation  of  the  inverse  problem  (Section  1),  we  now 
assume  that 


c(x)  »  1 


X  i  0. 


Lemma  A. 


Proof;  K,  being  the  smooth  part  of  a  solution  of  the  wave  equation 


must  solve  it  in  the  interior  of  the  light  cone.  For  x  O,  T(x) 


and  (3.7)  shows  that 


Finally 


A  solution  to  the  problem 


together  with  the  conditions  (3.8)  and  (3.9)  in  the  region  f-t  s  ^  ik  0 


However,  since  one  of  the  boundaries  of  this  region  is  characteristic,  the 
solution  is  unique  (see  [11],  Ch.  V).  q.e.d. 

Now  X~(s)  =  -s,  and  set  X(s)  =  X^(s).  Then  (2.6)  reads  for 


(X(s))K(X(s),t)  + 


dy  c  (y)K(y,  s)K(y, t) 


(X(s))K(X(s),t> 


dy  F (y  +  t )F (y  +  s )  + 


dy  c"  (y)K(y, s)K(y, t) 


Now  set 


H(s,t)  =  5  F(t+s)-r®  dT  F(t--)F(s--) 


We  have  proved 

Theorem  B .  The  smooth  part  K  of  the  solution  u  of  the  singular 


initial  value  problem  (2.1)  satisfies  for  0  ^  s  ^  t 
H(s,t)  =  I  c"^/^(X(s))K(X(s),t) 

+  dy  c"^(y)K{y,s)K(y,t)  (3.10a) 

*^0 

and  is  related  to  the  coefficient  c  ^  the  transport  equation 

K(x,T(x))  =  i  c^/^(x)('*  c^/^(c^/^)"  (3.10b) 

^  '"O 

We  shall  refer  to  the  system  of  integral  equation  (3.10)  as  the  G  L 
system,  partly  because  it  expresses  the  group  law  for  the  Cauchy  problem 
for  Cj.,  and  partly  because  it  is  related  to  an  equation  discovered  by 
Gel'fand  and  Levitan  and  derived  in  a  completely  different  way  in  their 
fundamental  paper  [13].  For  application  of  the  methods  developed  here  to 
the  problem  considered  in  [6],  and  a  discussion  of  various  aspects  of 
inverse  problems,  see  [6], [7]. 

The  converse  of  Theorem  B  will  be  proved  in  Section  7 . 


§4.  A  Priori  Estimates 

In  this  section  we  derive  some  necessary  conditions  for  scattering 
data.  We  shall  malce  use  of  the  substitution  x  s  =  T(x)  throughout.  W€ 
point  out  that  this  is  not  the  same  as  reduction  to  the  case  of  the 
Schrodinger  equation.  Indeed,  in  the  first  part  of  this  section,  leading 
up  to  the  estimate  (4.11),  the  Iternel  R  can  be  replaced  by  the  Iternel  J, 
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defined  in  the  r.ext  section.  The  operator  IK  can  be  replaced  by  an 
2  2-1 

operator  JI  :  L  (dt)  -•  L  (c  (x)dx),  and  all  mention  of  the  coordinate 
transformation  x  m  s  can  be  eliminated.  In  higher-dimensional  problems, 
the  volume  element  of  the  Riemannian  metric  associated  with  the  relevant 
hyperbolic  p.d.e.  will  play  the  role  of  c~^{x)dx,  so  our  methods  conform 
to  the  rubric  laid  down  in  Section  1,  Also,  the  last  part  of  the  section, 
leading  up  to  the  bounds  (4.16),  depends  only  on  x  «-•  s  as  arc-length 
parameterization  of  the  geodesics  of  the  above-mentioned  metric,  so  this 
is  again  an  admissible  trick. 

Note  that  the  kernel  H  is  symmetric.  It  follows  that  equation  (3.10a) 


also  holds  with  s  and  t 

Now  define 

interchanged  for 

0  i  t  ^  s. 

K(s,t)  =  2c"^/^(X(s))K(X(a),t)  for 

0  ^  s  i  t 

(4.1) 

K(s,t)  =  0 

for 

s  >  t  ^  0 

Since  c”^(y)dy  =  ds  with 

(3.10)  as 

s  =  T(y),  i.e.. 

y  =  X(s),  we 

can  re-write 

4H(s,t)  = 

K(s,t)  +  r  dx  K(t, 

•'o 

S)K(T,t) 

(4.2) 

2 

For  T  >  0  denote  by  IK  the  Volterra  operator  on  L  [0,T]  defined  by 


Kcp  (t  )  =  dt  K(T,t)ip(t) 

T 

2 

for  c  6  L  [0,T].  Denote  by  IH  the  symmetric  operator  with  kernel  H. 

The  hypotheses  on  c  are  sufficient  to  ensure  that  IH  is  a  Hilbert- 

2 

Schmidt  operator  on  L  [O, T].  Denote  finally  by  H  the  identity  operator 
2 

on  L  [0,T]  (with  kernel  {(r-t)).  Then  (4.2)  can  be  written 

n  +4  M  -  (n  +  IK  (n  +  K) 


r 


(see  [6],  eqn.  4.10,  also  [14],  eqn.  8.1,2).  This  shows  that  H  +  4  H 
must  te  positive  definite  (since  H  +  3K  is  invertible).  According  to 
the  Fredholm  character  of  H  +  4  M ,  we  must  in  fact  have 

n  +  4  m  ^  e(T)  >  o.  (4.3 


Obviously  t(T)  is  a  monotone  nonincreasing  function  of  T.  We  note 


for  later  use  the  identity 


€(T)“^  -  sup{'l(Il+  ]K)"^(cp)|| 


L‘'ro,T]  L''[0,T] 


which  follows  immediately  from  (4.3). 

Next,  we  extract  from  (4.3)  some  a  priori  estimates  which  will  be 
crucial  for  the  next  section.  Denote  by  the  functions 

H^.(s)  =  H(s,t),  Kj.(s)  »  K(s,t) 

and  note  that  (4.2  )  may  be  written 

4H^  (s)  =  (  (n  +  K  )'^  K^.)  (s),  O  ^  s  ^  t. 


It  follows  that 

16  H.  ^  16l'H.|l^,  =  <K.,  (n  +  K  )  (n  +  1K^)K^>  -  (4.5) 

l/[0,T]  ^  L^[0,t]  L^[0,t] 

Now  the  invertible  self-adjoint  operator  (11+  IK  )  (H  +  IK  )^  has  the  same 
spectrum  as  the  operator  (11+  IK  )^  (n  +  IK)  =  11+4  H,  hence  in  particular 
the  same  lower  bound.  Therefore  you  may  combine  (4.3)  and  (4.5)  to  obtain 

^  €(T)llK.|l^.  (4.6) 

‘  l2[o,t]  *=  l2[o,t] 


Now  ' may  be  estimated  in  the  following  truly  crude  fashion 


^2("F!'^2 

L^[0,T]  L‘^[0,2T] 


+  T^11P11‘*. 


L‘[0,T] 
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Also 


IJ®  dT  K(T,s)K(T,t)  i  «  l<Kg,K^.>|  2  ^  2 


L"[0,tj 


L^[0,s]  L^[0,t] 


Hence,  recalling  (4.6)  and  using  (4.7) 


j  dy  c”^ (y)K(y, s)K(y, t) I 


i  8€(T)’^(I'F!:2 

L  [0,2T] 


+  T^|lFll'*, 


L  [0,T] 

Combined  with  (3.10a),  this  yields  the  estimate 

|i  c"^/^(X(s))K(X(s),t)  I  ^  I'Fl'  +i'F"^, 


L"[0,2T]  l"[0,T] 


+  8€(T)~^[:'f1'^, 

L  [0,2T] 


L  [0,T] 


(4.8) 


(4.91 


valid  in  the  range  0  £  s  ^  t  £  T. 

Finally  we  give  a  priori  Isounds  for  c, 

in  terms  of  T  and  K*  =  sup  120”^^^^  (x  (t) )  K  (x  (t)  ,  t)  |  .  Note  that  the  latter 

0£t£T 

quantity  has  just  been  estimated  in  terms  of  F. 

According  to  (3.10b),  for  0  £  x  £  X(T), 


4c 


Set 


Then 


*^^^(x)K(x,T(x) )  -  /  0^/2  (^1/2)- 

0 

=  i  c' (X)  -  j  f  c“^(c')^ 
^  0 


■1/2 


g(x)  -  8  c  (x)  K  (x,T(x)) 


c(x)~^  g(x)  «  (loqc)'(x)  -  j  c~^ (x)  J  c~^(c')^ 
Since  c(0)  >  1,  obtain 

X  X 

log  c(x)  >  /  c  ^  g  ^  -4K*  /  c  ^  i  -4K*T 

"0  0 


i  .e. 


c  >  exp(-4K*T) 


(4.10) 


(4.11) 
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which  is  the  required  lower  bound  for  c. 

To  obtain  upper  bounds  for  c,  write 

j  g(X(s))  =  K(s)  =  K(s,s),  c?  =  •  =  A 

Then  (4.10)  may  be  rewritten 

2K{t)  =  (log  cp) '  (t)  -  ((logcp)’)^. 

It  follows  that  cp  obeys  the  equation 

Cp  5=  Qcp 

where  Q  =  2K.  (This  relation  is  also  part  of  the  Liouville  reduction 
which  figures  in  other  treatments  of  this  inverse  problem  —  see  e.g.  [2]). 


Note  that 


cp(0)  =  1,  b(0)  -  0,  K(0)  =  O. 


Hence  cp  is  the  solution  of 

cp(t)  =  1  +  ds(t  -  s)ip(s)Q(s)  =  1+2'**^  ds  cp(s)K(s) 

“'o  ‘0 

-  21'^  ds(t  -  s)ip(s)K(s) 

"o 


(4.12a) 


Also  (J)  is  the  solution  of 


ip(t)  =  2cp(t)K(t)  -2J*^  ds  5!(s)K(s) 


(4.12b) 


Define  ^  =  ip  -  2Kcp.  Then  the  Volterra  system  (4.12)  may  be  rewritten  as 


cp(t)  =  1+J^  da  2K(s)  (1 +2(s -t)K(s))cC!(s) 
■f  ds  2K(s)  (s  -  t)(,(s) 


(4.13a) 
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4(t)  =  -4r*^  ds(K(s))^(p(s)  -2r*^  ds  K(s)t(s) 


(4.13b) 


Set 


K**  =  max{4K*  +8T(K*)^,8(K*)^,4TK*,4K*]. 

Then  the  following  estimate  is  easily  derived  for  the  solution  of  (4,13): 

•  * 


llcpjl  ^  exp  TK 

L"[0,T] 


(4.14a) 


Also 


L*[0,T] 


(4.14b) 


In  view  of  the  definition  of  f,  this  entails 

^  (1 +2K*)exp  TK**  (4.15) 

L'LO^T] 

Now  i  ~  ^  (c^^)  =  c  ^  “  2  c'.  Thus,  combining  (4.15), 

(4.14a)  and  (4.11),  we  get  the  estimates,  valid  for  0  ^  x  ^  X(T), 


** 


and 


exp(-4K  T)  i  lc(x)|  ^  exp(2TK  ) 
lc'(x) 1  ^  2 (1  +2K*)exp  T(K**  + 4K* ) 


(4.16a) 


(4.16b) 


§5.  Solution  of  the  GL  system 

With  this  section  we  begin  the  solution  of  the  inverse  problem  as 
stated  in  Section  1.  The  first  step  is  to  show  that  the  GL  system  as 
presented  in  Theorem  B : 
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H(s,t)  =  i  (X(s)  )K(X(s),  t)  +  dy  c‘^  (y)K(y,s)K(y,t) 

‘  o 

O^s^t^T  (5.1a) 

K(x,T(x))  =  3  c^/^(x)'’’‘  (5.1b) 

■'0 

has  a  unique  solution  fK,cj,  where  c  is  defined  on  0  ^  x  ^  X(T), 
and  K  is  defined  in  {(x>t)  :0<^X;^X(T),  T(x)  ^t;^2T-T(x)}. 

Since  the  domains  on  which  the  solutions  are  defined  are  themselves 
defined  by  part  of  the  solution  (namely  c),  the  problem  has  something 
of  the  nature  of  a  free  boundary  problem. 

We  note  that  continuous  solutions  are  trivially  unique,  in  view  of 
the  Vol terra  character  of  (5.1  ). 

The  system  (5.1)  will  only  have  a  solution  as  described  when  H  has 
the  positivity  property 

n  +  4  IH  ^  6(T)  >0  (5.2) 

in  the  notation  of  the  last  section,  for  the  reasons  explained  there.  Our 
goal  is  to  show  that  this  necessary  condition  is  also  sufficient. 

First  introduce  the  function 

J(x, t)  =  2c"^/^(x)K(x,t) 

and  rewrite  (5.1)  as 

4H(s,t)  =  J(X(s), t)  +  dy  c~^(y)J(y,s)J(y,t)  (5.3a) 

‘O 

J(x,T(x))  =  i  r’'  c^/^(c^/^)"  (5.3b) 

2  ^0 


We  have  from  (5.3b),  (4.10) 
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Here  T(x)=  r**®  c"^  +  for  Xq^x^x^+Ax.  Similarly, 

■0  “Xq 

J„{x.t)  =  H(T„_^(x),t)  -f  dy  c;^i(y)Jn.i(y.Vl(’'>>'"n-l<^'^> 

’'o 

for  Xq  ^  X  ^  Xq  +  Ax 
t  ^  O. 

Now  suppose  6  >  0,  and  select  c*  with  0  <  c*  ^  ~ 

claim  that,  for  Ax  small  enough,  we  have  Cj^(x)  ^  c*,  ^  x  ^  x+Ax. 

In  fact,  suppose  that  this  is  so  for  c^,  k  =  0, l,...,n-l  (it  is 

obviously  true  for  n  =  O).  As  the  first  part  of  the  induction,  we 

estimate  J_  , ; 

n— 1 

’'o 

Now 

X  , 

H(s,t)  -  4[F(s  +t)  +  dT  F(s  -  T)P(t  -  t)]  -  r  °  <3y  c"-^  (y)  J(y,  s)  J(y,  t) 

*^0  '^0 

So  long  as  2(^0  ^  which  we  assume  for  the  moment,  we  have 


lH(s,t)|  ^  11F| 


L"[0,2T] 


+  IlFl 


L^[0,T] 


Ihe  second  summand  on  the  r.h.s.  is  estimated  by  (4.10): 

ir  °  dy  c"^(y)J(y,s)J(y,t)|  ^  32  €(T)"^[  llFjl^,  + 

•'o  L^[0,2T] 


l'‘[o,t] 


So 


lH(s,t)|  ^  411F||  +  (4  +32€(t)"^ 

L"[0,2T] 

+  326{T)"^  T^||F||S 

L‘‘[0,T] 


[0,2T] 


=  N(F,T,  e) 


(5.7) 
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Corabineci  with  (5.6)  this  yields 


J„.l  ,  ^  N(F,T,  .)  +  (X  -Xq)c;^  J^_2  I 


where  J means,  for  the  moment,  supf  J(x,  t)  :  ^  x  +  Ax,  t  ^  O], 

Now  suppose  that 


Jn.2  »  ^  ■t‘')N(F,T,  t), 


Then 


-1 


Jn-1  [1  +  (X  -  Xq)c/{1  +  J  ')N(F,T,  ■  )  ]N(F,T,  () 


Suppose  that  Ax  is  so  small  that 

-1 , 


Ax  c»  (1  +  »  ')N(F,T,  (.)  s,.  ' 
Then  we  hav'o  once  again  tliat 

^^1-1  cc  ^  ^  ‘  ')N(F,T,  c) 

To  complete  the  induction,  notice  from  (5.5.c)  that 

k(x,x^^)  _  c(x^)  c*  + 


(5.8) 


(5.9) 


Tlien 


so  that  (from  the  definition  of  c  ) 

n 


Cn(x)  _  c*  +  ■  -  .1  dy  (Vl  )  ^  c*  t  •  -4(1  +  ^  ' )  N  (F,  T,  ,. )  Ax 


so  wc  have  proved 

Lemma  l<  .  Suppose  ,  '  '  O,  and 


Ax  ^  minlc*[  (1  f  ■  •)N(F,T. f)  ', f 4 (1  +  >  ' )N (F, T,  ^ )  1'^ ‘ . 


(X)  1  c*,  n  0, 1, 2,  .  . . 

Xo  ^  X  Xp  +  ^x. 


;  J 


r- 
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Note  that  ‘  '  >  O  is  arbitrary  here,  whereas  f  must  be  chosen  so  that 
0  next  step  begins  with  the  assumption  of  a  Lipschitz 

bounii  on  F : 


'F{t)  -  F(s)  I  L  s  -  t 


(5. 101 


It  follows  that 

H(  ,t)-II(  .s)  ^  F('4t)-F(’+s)|+l’'  do  F(--o)(F(t-o)-F(s-ol)l 

■  0 

V  L' s  -  t  I  (1  +  Tl'F  )  L,'s-tl  (5.11) 


Set  J,  ,(y)  *:  J(y,s)  -J(y,t).  It  follows  from  (5.5a)  that 

S  ,  t 

J,  ,(y)  =  4[M(T(y),s)  -H(T(y),t)l  +  di  c"^  ( i  )  J  (  t  ,  y )  J  (,) 

S  ,  I  .  Q  S,  I 

Tliis  is  a  linear  Voltorra  equation,  whence  follows  the  estimate  (for 
O  y  ^  min  (  X  (  s)  ,  X  ( t  )  ) 


\s.l  L  |0,x  ^1' 

0 

exp[x^^c“^  J  I  (  s  -  t  ) 


L‘''fO,T] 


=  s  -  I 


(5. 12) 


Here  J  '  =  supf  J(x,t)  jO  ^  x  ^  t  ^  0}  is  estimated  by  (4.11)  iii 

terms  of  F,  so  is  estimated  in  terms  of  F,  c*,  T, 

Next  we  observe  Dial 

'  o  o 

So 


'^'n"'’  -■''n-l*'''  i.  '^n-l‘y'  -^n‘^'’ I 

.  -2 

^  Axe*  c„  1  -  c„ 

*  n-1  n  *' 


(5.11) 


i  •  - 

MHMiaiiMaMlMyiiyttAiliK 
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where  for  the  movement 

i!c„  =*  S«P  ,(y)-c  {y)l 

and  so  on. 

According  to  the  definition  of  H, 

~  X  . 

H(s,t)  -  H(t,T)  -  4[H{s,  t)  -H(t,T)]  -  r  °  dy  c"'^  (y )  [  J  (y.  s  )  -  J(y,  t)]J(y, 

‘0 

Hence 


-H(t,T)|  £  4Lj1s  -  t|  +XpC;^|!j!'^  Ljls  -  t 


“  Ljls-tl 


(5.14) 


where  again  is  estimated  in  terms  of  F,  c*,  T,  and  x^.  Next 

estimate  for  ^  y  ^  x 

lJn<y,t)  -  Jj^(y,s)|  ^  |H(T^.l(y),t)  -H(T^_j(y),s) 

^  i/x^ 

i  Ljit  -si 

+  (y-Xo)c;^(!Jn_ii;^  sup  ~'^n-i 


(z, s)  J 
(5.15) 


According  to  (5.8),  (5.9), 


-1, 


i  6'. 


Select  so  that 

1-4  i  (1  -6')"^  Lj 

(which  introduces  the  new  restriction  6'  <  1).  Then  (5.15)  and  obvious 
induction  guarantees  that 


Liw.-. 


Note  that  c, 


We  suppose  that  d  >0  such  that 


+  Ax ,  where 


(d  )'‘x  +4(1  +  6)N(F,T,  6)  +  t 


Wien 


(d  )‘‘x  +4(1  +  6)N(F,T,€)  + 


where  we  have  used  (5.9)  and  Lemma  a.  So  we  have  proved 


Lemma  C.  Provided  Ax  satisfies  the  bounds  of  Lemma  B  and  addition 


we  have  the  estimates  for  all 


We  are  now  ready  for  the  main  estimates.  First 


(x))J 


(x))J 


3 
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Ax  Lj  c*^!'Cn-l  ■ 

-1  -  .  -1 


*  r  ‘^y  ‘=n-2‘y>iJn-i<y»Vi(’‘>  -  Vi‘y'V2<’‘’>MJn-i‘’ 


^  r  ‘‘y  C2<yHJn-l‘y'V2<’‘>>  -'^n-2‘y'V2‘’'>>l!'^n-l 


-1 


*  r  ‘‘y  <-2‘yMV2‘y'V2<’‘>viJn-i‘y'‘’  -  V2<y'‘’ 

’‘o 


Ax  c*^[L3  +  llJ„_j"^]liCj^_j  -  Cj^_2| 


+  Ax  c 


*  .•^"V2»JK-i- V2IL 


^  Ax  A"c  I  -  c  ,  ♦Ax  B''J„  1  ”  o 
^  n-1  n-2  •  n-1  n-2 


(5 


where 


A  -  c;^[Lj  ♦  (1  ♦  f  ')N(F,T,  ,)[  (1  ♦  ?')N(F,T,  €3  +  0;^  1  ] 

B  -  2c;^  (1  ♦  ‘  ')N(F,T,  O . 


Next  estimate 


s‘*’ ^  V  “y  Vi‘y'Vi‘y”  -  V2<y>  V2<y’>i 


*0  “o 


^  V  ‘'y  Vi‘y-Vi‘y’' -  Vi<y'V2<y'M 


♦  4  ciy  Vi(y.V2(y)  -  J„.2(y.V2(y))l 
*0 

"5.  “  ‘^y  .■'’  S-l‘^n-2l‘=n-l-S-2i 

O  ’‘0 


r-,  t)  I 
(y,t)| 


.18) 


U^mtkjlLrn 


sup  I Jn(y, s)  -  J  (y, t) I  ^  L. ! s  -  t ' ,  n=  0,1,2,...  (5.16 

^O^S^o+Ax 


Note  that  c'  satisfies 
n 


’‘o 

Hence, 

°  ’‘o 

* 

We  suppose  that  d  >0  such  that 


l"[o,Xq] 

** 

and  lc^(x)  j  <,  d  ,  k  =  0,1,  ..  .,n  -  1,  ^  x  ^  +Ax,  where 

d**  =  i  c;^(d*)^XQ  +4(1  +  6)N(F,T,  €)  +  6. 


|Cn(x)!  ^  5  c;^(d*)^XQ +4(1  +  6)N(F,T,€)  +  |  Ax  c;^(d**)^ 

where  we  have  used  (5.9)  and  Lemma  A.  So  we  have  proved 

Lemma  C.  Provided  Ax  satisfies  the  bounds  of  Lemma  B  and  addition 


Ax  i  2c*  (d  )  6 

we  have  the  estimates  for  all  n  ^  O: 

|Cn(x)l  i  d**,  Xq  i  X  i  Xq  +AX 
we  are  now  ready  for  the  main  estimates.  First, 

^  If  ‘‘y(^;-i<y)Vi<y'Vi< 

’‘o 


^n-2  ^y^'^n-2  ^ 
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1  pX 


dy 


-1 

'n-2 


"n-l 


1  4(Ax)2  c;2  L4llc„.i-c„.2lL+4Axl<J^_j 
+  J(Ax)2  c:2(d**)2|,c^_^_c^_^||^ 

+  ^(Ax)2  c:^ 


"n-2 1 


'n-2 


O' 


(5.19) 


Finally, 

l'=n<^'  -‘=n-l<’'H  ^  41j^_^(x,T„_i(x))  -J^.2(x,T„.2(x))| 

!‘=;-l<‘=n-l^'-‘=;'2<-n-2)'l 

’‘o 

^  4|J„_i(x,Vl(x))  -J„.;^(x,T„_2(x))| 

'‘!’^n-l^’‘»'^n-2^*^  ^  ■ ‘^n-2^’‘*’’n-2  I 

5  J  ^n-l‘^n-2l‘^n-l  * ‘^n-2l  ^*^11-1^ 

’‘o 

*  k  ‘^n-2  I  ‘^n-1  *  ^^11-2  !  I  ‘^n-l  “  ‘^n-2  I 

i  4AX  c;2|!c^_^-c„_2ll^^4l!J„_^-J„_2l|^ 

+  5  Ax  c*  (d  )  l|Cn_i  - 

+  Ax  c;^  ‘***r«'n-l  ■‘'n-2lL  <^*20) 


The  estimates  (5.18),  (5.19)  and  (5.20),  ta)ten  together,  show  that, 
provided  that  Ax  >  0  satisfies  the  bounds  in  Lemmas  b  and  c>  and  is 
possibly  smaller  yet,  the  sequences  {J^)/  f^^n^  converge 

uniformly  on  x^  ^  x  ^  Xq  +Ax  to  solutions  of  (5.5).  Ihe  numbers  which 
determine  how  small  Ax  must  be  are  c*,  f ',  N(F,T, €),  j,  L  (in  (5.10)), 
F  ,  T,  "jll  ,  d*,  and  X-.  Of  these,  N(F,T,  6),  L,  I'Fl’  and  T  are 

CP  er  ^  '  flB 

determined  by  the  data,  c*  is  estimated  from  below  by  (4,16a),  d  is 


I 


1 


f 
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^  Ax  Lj 

f  <ay!<=n-l(y) 

*0 

’‘o 

+  f  =;;-2<yHJn-l<i''V2(*>^  -Jn-2<y'V2‘*)>llJn-l<y'*=n 

’‘o 

^  f  '"y  <=;-2<>'>!'^n-2(y'V2(^))MJn-l<i''">-'"n-2(y»">l 

’'o 

^Axc;2[L3  +  llJn.il'^]K.l-c„.2lL 

+  AX  c:^||J^_3!1^  ^4K-l-‘=n-2lL 

^  AX  c:^(l!Jn-llL-^K-2llJl|Jn-l- V2L 

^  Ax  Al!c^_j^-c^_2lL+Ax  b|!j„_3- J^_2''» 

where 

A  =  c;^[L3  +  (1  +  6')N(F,T,  e)[  (1  +  6')N(F,T,  €)  +  cZ^  L^]] 

B  =  2c;^(l  +  6')N(F,T,  6). 


Next  estimate 


i‘=n<’‘’ ^  -'^n-2<y*V2<y>n 


0 


-1 


2  _-l 


O 

^  '*f  ‘'y(Jn-i(y*V2(y>  - Jn-2(y'V2(y)’i 


*  2  J  ®n-l‘^n-2l®n-l  ■  ^n-2l 

0  ’‘o 


■kuikdiflii 


jm.-jr.i^  .-  .I,.. 


■  ■.■AV.KJt.-:-.  1..  .  Jr... 


^  4  (Ax)  c 


Finally 


The  estimates  (5.18),  (5.19)  and  (5.20),  talcen  together,  show  that 


provided  that  Ax  >  O  satisfies  the  bounds  in  Lemma s  b  and  c»  and  is 


possibly  smaller  yet,  the  sequences 


x_  +Ax  to  solutions  of  (5.5).  The  numbers  which 


uniformly  on  x 


determine  how  small  Ax  must  be  are  c 


determined  by  the  data 


c^,  is  estimated  from  below  by  (4.16a) 
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estimated  by  (4.16b),  I'J  ^  is  governed  by  the  main  a-priori  estimate 
(4.11),  and 

Xq  i,  c*T 

* 

where  c  is  an  upper  bound  for  c,  given  by  (4.16a).  It  follows  that, 

for  given  T  and  F  satisfying  the  positivity  condition  (5.2),  Ax  may 
be  chosen  independently  of  so  long  as  x^  X(T).  Thus  finitely 

many  repetitions  of  the  iteration  scheme  outlined  above  suffice  to  deter¬ 
mine  c  on  the  interval  [0,X(T)],  and  J  on  the  corresponding  domain. 

The  system  (5.5),  (and  with  it  (5.1))  has  therefore  been  solved,  as 
promised. 

By  differentiating  the  G-L  system  (5.1),  one  obtains  systems  of 
Volterra  equations  for  the  derivatives  of  c  and  the  partial  derivatives 
of  K.  These  are  (essentially)  linear  systems.  Without  carrying  out 
the  details,  we  state  that  these  systems  for  the  derivatives  possess 
continuous  solutions.  As  in  the  case  of  the  GL  equation  itself,  c  winds 
up  with  one  more  derivative  than  K,  and  K  has  as  many  derivatives  as  F. 
The  solution  of  the  GL  system  therefore  defines  a  map  F  {K,c).  It  is 
easily  verified  that  the  positivity  condition  is  stable  under  perturbation. 
It  follows  that  the  map  F  i-*  {K, c}  is  continuous  in  the  obvious  sense  of 
c"’  norms,  as  outlined  in  the  introduction, 

§6.  The  lower  bound  on  e 

For  reasons  explained  in  the  introduction,  we  now  estimate  €(T) 

*  *  * 

(see  (4.3))  in  terms  of  the  bounds  c„,c  ,d  and  e  which  we  suppose 
given: 


f. 
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c*  1  c(x)  ^  c  , 

)c'(x)  I  ^  d*, 

|c''(x)  ^  e*,  O  ^  X  ^  X(t). 

Recall  that  K  selves  the  boundary  value  problem 

2  2 

(x)  -2^\k(x,  t)  =  O,  t  ^  T(x) 
hx  I 

K(x,T(x))  =  i 

=  g  c^/^(x)c'(x)  -  c^/^(x)  c~^(c')^ 

(X  N  O) 

Also,  according  to  Lemma  C  (Sec.  3)  (recall  c  =  1  for  x  <  O) : 


(6.1) 


K(x,-x)  =0,  X  <  O  . 

It  follows,  as  in  [11],  Ch.  V,  that  K  is  the  solution  of  an  integral 
equation  of  Volterra  type.  In  fact,  if  you  denote  by  C(x, t)  the 
intersection  of  the  backward  light  cone  with  vertex  (x, t)  with  the 
forward  light  cone  with  vertex  (0,0),  and  by  r(x, t)  the  x-coordinate 
of  the  intersection  of  the  characteristic  curve  of  negative  shape  through 
(x, t)  with  the  characteristic  curve  of  positive  slope  through  (0,0), 
you  eventually  obtain 

5  '■  *  5 

+  |[c'(r(x,t) 

(recalling  that  J(x, t)  =  2c^/^ (x)K(x,  t) ) .  Now 

vol  C(x, t)  ^2  c*t^. 

It  follows  easily  that 

I  J(x,t)  I  ^  Pf  (  ^  4*Q  t) 


(6.2) 


7 
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where  f  is  the  entire  function  with  power  series 

8 


f(2)  =  +  ^  .  5-3 


4  6 

z  ^  ^ 


7.5.3 


.} 


and 


P  =  sup  i[c'(x)  -  i  c"^(c')^] 
O^Xi;X(T)  “  ‘0 

1  *  1  ♦  -  1  *  ? 

^  i[d  +i  c  c/(d  )^T] 

Q=  sup  [i  c*  (X)  +3  c"^(x)  (c'(x)  )2] 
O^x^X(T) 

*>2 


^  5  ®  +2  c*  (a  ) 


2  4 

Recalling  the  definition  (4.1)  of  K  and  of  the  operator  IK,  one  sees 

that  K  obeys  the  same  sup  norm  bound  as  J,  namely  (6.2)  above.  The 

-1  ^ 

Iternel  of  the  inverse  operator  (H  +  K )  ,  say  K,  is  then  easily  bounded: 


"K'l„  i|KlL  exp(;|K!'  T) 


(6.3) 


It  follows  immediately  that 

|;(ll+  lK)"^cs!!^  ^  (1  +  !1k!'^t)^||cp|,^ 

So  (see  (4.4)) 

e(T)~^  ^  (l  +  l|K’i^T  exp(llKll^T)  )^ 

which  together  with  the  bound  (6.2)  estimates  6(T)  in  terms  of  T  and 
the  a  priori  information  on  c,  as  desired. 


§7 .  Equivalence  of  GL  and  Chudov  systems 

We  show  by  direct  computation  that  the  solution  of  the  GL  system 
constructed  in  Section  5  solves  the  Chudov  boundary  value  problem 

=  0  in  r t  ^  0,  0  ^  X  ^  X(t)}  (7.1a) 

'  dt  3X  / 

K(0,t)  =  F(t) 


Now  it  was  shown  in  §2  that  the  smooth  part  of  the  solution  of  the  initial 
value  problem  (2.1)  solves  (7.1).  Since  the  solution  of  (7.1)  can  easily 


be  shown  to  be  unique,  it  follows  that  the  solution  of  the  GL  system  is 
in  fact  the  smooth  part  of  the  solution  of  (2.1)  with  corresponding 


hence  that  we  have  solved  the  inverse  problem 


First  note  from  (3.10c)  that 


the  solution  of  the  GL  system  as  constructed 


in  fection  5.  We  assume  of  course  that  F  satisfies  the  positivity 


condition  (5.2)  and  is  of  class  C  on  its  interval  of  definition  so  that 


K  and  c  are  of  classes  C  and  C  respectively  on  their  domains  of 


definition 


Setting  s  =  O  in  (5.10)  and  recalling  the  definition  (3.lOc)  of 


Using  the  definition  (3.l0c)  again,  and  the  requirement  F(0)  =0,  one 


sees  that 


On  the  other  hand,  from  (5.1)  one  obtains  (recalling  c(0)  *  1,  c'(0)  =  O 


I 


Thus  K  obeys  the  boundary  conditions  on  the  Chudov  system.  It  romuins 


only  to  verify  that  K  solves  the  wave  equation  in  the  interior  of 


the  light  cone.  To  do  this  first  compute 


(X(s)  )D,K(.X(s),  t) 


dy  t  (y)K(y, s)D,K (y, t )  ] 


(X(S))  (c'(X(s)  )  )''k(X(s),  t) 


(X(s) )c" (X(s) )K(X(s),  t) 


(X(s) )D,K(X(s),  t) 


[c“'^  (X(s)  )K(X(s),  s)K(X(s),  t)  ] 


c"-"  (X(s)  )D  K(X(s),  s)K(X(s),  t) 


Using  (7.1c)  one  sees  that 


(X(s) )K(X(s),s) 1 


(X(s)) 


K(X(s),s)  -  c  ■^(X(s))c'(X(s))K(X(s),s) 


(X(s)l (c'(X(s))) 


(X(s)  )c''  (X(s)  ) 


tc~'^(X(s)  )K(X(s),s)K(X(s),t)  ]  +  c"  (X(s))D-K(X(s),s)K(X(s),t) 


c"-^(X(s)  )c'(X(s)  )K(X(s),s)K(X(s),t) 


+  2c“-"(X(s)) 


K(X(s),s)K(X(s),t)  -  D,K(X(s),s(K(X(s),t) 


+  K(X(s),s)Dj^K(X(s),  t) 


(X(s))  (c'(X(s1))‘-lK(X(s),t) 


(X(s))c  (X(s)) 


Dj^K(X(s),  s)K(X(s),  t)  +  K(X(s),  s)DjK(X(s),  1) 
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Hence 


- ^jH(s,t)  -  I  c"^/^(X(s))[D2K(X(s),t)  -  c^(X(s))D^K(X(s),t) 

+  dy  c”^(y)[K(y,s)D2K(y,  t)  -D2K(y,s)K(y,t) 

'0 

+  Dj^K(X(s),  s)K(X(s),t)  -  K(X(s),s)D^K(X(s),t)  (7.31 

Now  add  to  (7.3)  the  identity 

0  =  -D^K(X{s),s)K(X(s),t)  +  K(X(s),s)Dj^K(X(s),t)  +  D^K  (0,  s  )K  (0,  t ) 

-  K(0,  s)Dj^K(0,  t) 

-  dy  c“^  (y ) [K(y, s) (C2 (y )DjK(y, t ) ) 

•  0 

-  (c^(y)D^K(y,s))K(y,t)  ] 

obtained  by  integration  by  parts.  After  a  little  manipulation,  one  gets 

- ^)H(s,t)  =  D,K(0,s)K(0,t)  +  K(0,s)D,K(0,  t) 

-=  2  (X(s)  )  [D2K(X(s),t)  -  (X(s)  )D^K(X(s),  tl  ' 

+  dy  c~^(i):K(y,9){D^K(y,t)  -  (y  )D?K  (y,  1 11 

-  {D2K(y,s)  -  c^(y)D^K(y,s))K(y,t)]  (7.41 

Comparing  with  (7.2)  and  using  (7.1b),  one  sees  that  the  l.h.s.  of  (7.41 

in  fact  vanishes  identically.  Eq.  (7.4)  is  therefore  a  linear  integral 

equation  of  Volterra  type,  since  K  is  continuous  and  c  >  O,  for 
2  2  2 

D2K  -  c  Dj^K.  It  follows  that  this  expression  vanishes  in  the  domain  covered 
by  the  limits  of  integration,  that  is,  in  the  light  cone.  This  result 
shows  that  the  solution  K  of  the  GL  equation  solves  the  Chudov  system 
also.  Since  we  showed  in  Sections  2  and  3  that  the  solution  of  t)ie  Chudov 
system  solves  the  GL  system,  it  follows  that  these  two  sy-. terns  are 


J8 


completely  equivalent.  It  also  follows  that  c,  as  part  of  the  solution 
of  the  GL  system  for  prescribed  F,  solves  the  inverse  problem  stated 
in  Section  1,  Together  with  the  stability  statement  of  Section  5,  this 
constitutes  a  proof  of  Theorem  A . 
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